STRONG SOLUTIONS TO THE THREE-DIMENSIONAL 
COMPRESSIBLE VISCOELASTIC FLUIDS 



XIANPENG HU AND DEHUA WANG 

Abstract. The existence and uniqueness of the local strong solution to the three- 
dimensional compressible viscoelastic fluids near the equilibrium is established. In addi- 
tion to the uniform estimates on the velocity, some essential uniform estimates on the 
density and the deformation gradient are also obtained. 



1. Introduction 

Viscoelastic fluids exhibit a combination of both fluid and sohd characteristics, and 
keep memory of their past deformations. The interaction between the microscopic elastic 
properties and the macroscopic fluid motions leads to the rich and complicated rheological 
phenomena in viscoelastic fluids, and also causes formidable analytic and numerical chal- 
lenges in mathematical analysis. We consider the following equations of three-dimensional 
compressible flow of viscoelastic fluids O [lOl [16] : 

pt + dw{pu) = 0, (1.1a) 

{pu)t +div{pu(g)u) - pAu - {X + p)Vdwu + VP{p) = div(/3FF'^), (1.1b) 

Ft-Fu- VF = VuF, (1.1c) 

where p stands for the density, u G the velocity, and F S M^^^ (the set of 3 x 3 matrices) 
the deformation gradient. The viscosity coefficients p, A are two constants satisfying p > 
0,2p + 3\ > 0, which ensures that the operator —pAu — (A + /i)Vdivu is a strongly elliptic 
operator. The pressure term P{p) is an increasing and convex function of p for p > 0, 
and in particular, we take P{p) = p^ with 7 > 1 a conatant. The symbol ® denotes the 
Kronecker tensor product, F^ means the transpose matrix of F, and the notation u • VF 
is understood to be (u • V)F. As usual, we call equation (|l.la|) the continuity equation. 
For system p.ip . the corresponding elastic energy is chosen to be the special form of the 
Hookean linear elasticity: 

W{7) = \\Y\\ 

which, however, does not reduce the essential difficulties for analysis. The methods and 
results of this paper can be applied to more general cases. 

In this paper, we consider equations (jl.ip subject to the initial condition: 

(p,u,F)|f=o = (po(a:),uo(x),Fo(x)), xGM^ (1.2) 
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and we are interested in the existence and uniqueness of strong solution to the initial- value 
problem (|l.ip - (jl.2p near its equilibrium state in the three dimensional space M^. Here the 
equilibrium state of the system (ll.ip is defined as: p is a positive constant (for simplicity, 
p = 1), u = 0, and F = / (the identity matrix in M^^^). We introduce a new unknown 
variable E by setting 

F = I + E. 

Then, (jl.ip becomes 

pt + div(pu) = 0, (1.3a) 
{pu)t + div (pu (g)u)- pAu -{p + A)Vdivu + VP{p) = dw{p{I + E){I + E^), (1.3b) 
Et + u-VE = VuE + Vu, (1.3c) 

with the initial data 

{p,u,E)\t=o = {po{x),uo{x),Eo{x)), xeR\ (1.4) 

By a strong solution to ()1.3p - (|1.4p . we mean a triplet {p,u,E) satisfying (jl.3p almost 
everywhere with initial condition (jl.4p . in particular, u{-,t) £ W"^'"^ and {p{t, ■), E{t, ■)) G 
W^''^ with q G (3,c«) in some time interval [0,T] for T > in this paper. 

When the density p is a constant, system (jl.ip governs the homogeneous incompressible 
viscoelastic fluids, and there exist rich results in the literature for the global existence of 
classical solutions (namely in or other functional spaces with much higher regularity); 
see O [6l [HI [131 [El [El [13 [13 [21] and the references therein. When the density p is not a 
constant, the question related to existence becomes much more complicated and not much 
has been done. In [15] the authors considered the global existence of classical solutions 
in of small perturbation near its equilibrium for the compressible viscoelastic fluids 
without the pressure term. One of the main difficulties in proving the global existence 
is the lacking of the dissipative estimate for the deformation gradient and the gradient 
of the density. To overcome this difficulty, for incompressible cases, the authors in |14j 
introduced an auxiliary function to obtain the dissipative estimate, while the authors in 
[I6j directly deal with the quantities such as Au + divF. Those methods can provide them 
with some good estimates, partly because of their high regularity of (u, F). However, in 
this paper, we deal with the strong solution with much less regularity in W^'"^, q G (3, oo], 
hence those methods do not apply. For that purpose, we need a new method to overcome 
this obstacle, and we find that a combination between the velocity and the convolution 
of the divergence of the deformation gradient with the fundamental solution of Laplace 
operator will develop some good dissipative estimates which may be very useful for the 
global existence. The local existence is established using a fixed point theorem. Some 
uniform estimates on the solution are also obtained. These estimates are essential for the 
global existence although one of the estimates needs to be improved in order to establish 
the global existence. 

The viscoelastic fluid system (II. ip can be regarded as a combination of compressible 
Navier-Stokes equations with the source term div(pFF^) and the equation (jl.lcp . For 
the global existence of classical solutions with small perturbation near an equilibrium for 
the compressible Navier-Stokes equations, we refer the reader to |22l [23l |24l |26] and the 
references cited therein. We remark that, for the nonlinear compressible inviscid elastic 
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systems, the existence of solutions was established by Sideris-Thomases in [29] under the 
null condition; see also [27] for a related discussion. 

The existence of global weak solutions with large initial data of (II. Ih is still an out- 
standing open question. In this direction for the homogeneous incompressible viscoelastic 
fluids, when the contribution of the strain rate (symmetric part of Vu) in the constitutive 
equation is neglected, Lions-Masmoudi in [20] proved the global existence of weak solu- 
tions with large initial data for the Oldroyd model. Also Lin-Liu-Zhang in [17j proved the 
existence of global weak solutions with large initial data for the incompressible viscoelastic 
fluids when the velocity satisfies the Lipschitz condition. When dealing with the global 
existence of weak solutions of the viscoelastic fluid system (jl.ip with large data, the rapid 
oscillation of the density and the non-compatibility between the quadratic form and the 
weak convergence are two of the major difficulties. 

The rest of the paper is organized as follows. In Section 2, we recall briefly the com- 
pressible viscoelastic fluids from some basic mechanics and conservation laws. In Section 

3, we state our main results, including the local existence and uniqueness of the strong 
solution to the system (jl.3p - ()1.4p . as well as uniform estimates which may be very useful 
for the proof of the global existence. In Section 4, we prove the local existence via a fixed- 
point theorem. In Section 5, we prove the uniqueness of the solution obtained in Section 

4. In Section 6, we establish some uniform a priori estimates, especially on the dissipation 
of the deformation gradient and gradient of the density. In Section 7, we improve the 
uniform estimates in Section 6. 



2. Background of Mechanics for Viscoelastic Fluids 

To provide a better understanding of system (jl.ip . we recall briefly some background 
of viscoelastic fluids from mechanics in this section. 

First, we discuss the deformation gradient F. The dynamics of a velocity field u{x,t) 
in mechanics can be described by the flow map or particle trajectory x{t,X), which is a 
time dependent family of orientation preserving diffeomorphisms defined by: 

^x{t, X) = u{t, x{t, X)), x(0, X) = X, (2.1) 
at 

where the material point X (Lagrangian coordinate) is deformed to the spatial position 
x{t,X), the reference (Eulerian) coordinate at time t. The deformation gradient F is 
defined as 

Ht,x) = —{t,x), 

which describes the change of configuration, amplification or pattern during the dynamical 
process, and satisfies the following equation by changing the order of differentiation: 

dHt,X) _ du{t,xit,X)) 

dt ~ dx ■ ^ ' 

In the Eulerian coordinate, the corresponding deformation gradient F(t, x) is defined as 

F{t,x{t,X))=F{t,X). 
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Equation (j2.2p . combined with the chain rule and (j2.ip . gives 

dtHt, xit, X)) + u • VF(t, x{t, X)) = dtF{t, x{t, X)) + ^^^^'f^'^^^ ■ ^^1^ 

dF{t,X) dvi{t,x{t,X)) du{t,x{t,X)) dx 



dt dX dx dX 

du{t,x{t,X)) ~ 
= F(t,X)=Vu.F, 

which is exactly equation (jl.lcp . Here, and in what follows, we use the conventional 
notations: 

du- dF- 
= 9^' = (^'^^^''^kj, (u • VF),j =UkQ^, 

and summation over repeated indices will always be well understood. In viscoelastic fluids, 
(jl.lcp can also be interpreted as the consistency of the flow maps generated by the velocity 
field u and the deformation gradient F. 

The difference between fiuids and solids lies in the fact that, in fiuids, such as Navier- 
Stokes equations |2l] , the internal energy can be determined solely by the determinant part 
of F (equivalently the density p, and hence, (jl.lcp can be disregarded); while in elasticity, 
the energy depends on all information of F. 

In the continuum physics, if we assume that the material is homogeneous, then the 
conservation laws of mass and of momentum become [3 [T^ [25| [27] : 

dtp + div(pu) = 0, (2.3) 

and 

dtipu) + div(pu u) - pAu - {p + A)Vdivu + VP(p) = div((det F)"^5f'^), (2.4) 
where 

pdetF = l, (2.5) 

and 

dW 

Sm = ^. (2.6) 

Here S, pSF~^ , W{F) denote Piola-Kirchhoff stress, Cauchy stress, and the elastic energy 
of the material, respectively. Recall that the condition (j2.6p implies that the material is 
hyperelastic pi] . In the case of Hookean (linear) elasticity [HI [HI [18] , 

W{F) = i|F|2 = ^tr(FFT), (2.7) 

where the notation "tr" stands for the trace operator of a matrix, and hence, 

S{F) = F. (2.8) 

Combining equations (j2.ip - (j2.8p together, we obtain system (jl.ip . 
If the viscoelastic system (jl.ip satisfies 

div(poF([) = 0, 

initially at t = (with Fq = / + Eq), it was verified in [16] (see Proposition 3.1) that this 
condition will insist in time, that is, 

div(/9(t)F(t)^) = 0, for t > 0. (2.9) 
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Another hidden, but important, property of the viscoelastic fluids system (jl.ip is con- 
cerned with the curl of the deformation gradient (for the incompressible case, see |131ll4j). 
Actually, the following lemma says that the curl of the deformation gradient is of higher 
order. 

Lemma 2.1. Assume that (jl.lcp is satisfied and (u,F) is the solution of the system (jl.ip . 
Then the following identity 

FikViFij = FijViFik (2.10) 
holds for all time t > if it initially satisfies (j2.10p . 

Again, throughout this paper, the standard summation notation over the repeated index 
is always adopted. 

Proof. First, we establish the evolution equation for the equality FikViFij — FijViFik- In- 
deed, by the equation (jl.lcp . we can get 

dNiFij + u • VV/Fjj + Viu ■ VFij = VmUiV/F^j + V/VmUiFmj- 

Thus, 

FikidtViFij + u • VViFij) + F/fcV/u • VF^ = FikVmUiViFmj + FikViVm'H-iFmj- (2.11) 
Also, from (jl.lcp . we obtain 

ViFijidtFik + u • VFik) = ViFijVmniFmk- (2.12) 
Now, adding (f2ll]) and (|2J2]) . we deduce that 

dt{Fik^lFij) + u • V(FifcV;Fij) = -FikViVi ■ VFij + FikVmUiViFmj 

m'^iFnij + ^iFijVm'^lFmk (2-13) 
= FikVm'^iViFmj + FifcV/VmUjFmj. 

Here, we used the identity which is derived by interchanging the roles of indices / and m: 

FzfcV/U • VFij = FikViUjnVmFij = ViFijVmUiFmk- 

Similarly, one has 

dt{FijViFik) + u • V{FijViFik) = FijVm^iViFmk + FijViVm^iFmk- (2.14) 
Subtracting (l2H]l from (j2J3]l yields 

dt{FikViFij - FijViFik) + u • V{FikViFij - F/jV/Fjfc) 

= VmUi(F;fcViFmj - F/jV/Fmfe) + VmUj (F^jF/fc - F^feF/j). 

Due to the fact 

VlVmUi = VmV/Uj 

in the sense of distributions, we have, again by interchanging the roles of indices / and m, 

VlVmUi{Fm.jFlk — FjnkFlj) = V iV m'^iFmjFlk — Vm,UjFmfcF/j 

= V/VmUjFmjFjfc — VmViUjF^fcFmj 
= (ViVmUi - VmVlUi)FikFmj = 0. 

Prom this identity, equation (j2.15p can be simplified as 

dt{FikViFij - FijViFik) + u • V{FikViFij - FijViFik) 

= VmUi(F;fcViFmj - FijViFmk)- 
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Multiplying (I2J6I1 by FikViF^ - F^jV/Fifc, we get 

dt\Fik^lFij - F/jVzFjfcp + u • \7\FikViFij - F/jV/Fjfc|^ 

= 2{FikViFij - FijViFik)VrrM^lk^l'Pmj - FijViF^k) (2.17) 

< 2||Vu||ioo(K3)A^^ 
where M. is defined as 

M = m&x{\F ikViFij - FijViFikl"^}. 

Hence, (|2.17p implies 

dtM + u-VM <2\\Vu\\Lo.(^^i^M. (2.18) 
On the other hand, the characteristics of dtf + u • V/ = is given by 

^X{s) = u{s,X{s)), X{t) = x. 
as 

Hence, (I2.17P can be rewritten as 

^ < B{t, y)U, U{0, y) = Mo{y), (2.19) 

where 

U{t,y) = M{t,X{t,x)), B{t,y) = 2||Vu||i^(K3)(t, X(t, y)). 
The differential inequality (j2.19p implies that 

ft 



Hence, 



Uit,y)<U{0)exp(^J^ B{s,y)ds^ . 
M{t,x) < A^(0)exp (^J 2||Vu||ioo(R3)(s)ds^ . 



Hence, if A^(0) = 0, then A4{t) = for all t > 0, and the proof of the lemma is complete. 

□ 

Remark 2.1. Lemma l2.ll can be interpreted from the physical viewpoint as follows: for- 
mally, the fact that the Lagrangian derivatives commute and the definition of the defor- 
mation gradient imply 

~ dXkdXj ~ dXjdXk ~ 

which is equivalent to, in the Eulerian coordinates, 

FikViF,j{t, x{t, X)) = FijViFik{t, x{t, X)), 

that is, 

F«fcV;Fjj(t,x) = FijViFik{t,x). 

Using F = I + E, (j2J0|) means 

"^kEij + ElkViEij = VjEik + EijViEik- (2.20) 

According to (12.20p , it is natural to assume that the initial condition of E in the viscoelastic 
fluids system (|1.3p should satisfy the compatibility condition 

VkEiO)ij + E{0)ikViE{0)i, = VjE{0)ik + EiO)ijViEiO)ik. (2.21) 
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Finally, if the density p is a constant, then we have the following equations of incom- 
pressible viscoelastic fluids (see [H [131 El HSl [13 [H] and references therein) : 

divu = 0, 

5tu + u- Vu-pAu + VP = div(FF^), (2.22) 
^tF + u- VF = VuF. 

For more discussions on viscoelastic fluids and related models, see [H [5l [71 [H [lOl [12l [171 
[201 [211 [23 [30] and the references cited therein. 

3. Main Results 

In this Section, we state our main results. The standard notations for Sobolev spaces 
W^''^ and Besov spaces Bp^ ([3]) will be used. Throughout this paper, the real interpolation 
method (0) will be adopted and the following interpolation spaces will be needed 

(L5(M3),t^2,,(j^3))^_ ^2(1-1)^ (L«(m3),W/1.«(m3))^_ bI;\ p,q>l. 

Now we introduce the following functional spaces to which the solution and initial con- 
ditions of the system (11. 3p will belong. Given 1 < p, g < oo and T > 0, we set 
Qt = R^ -K (0,r), and 

H;P''?(0,r) := {u : u G VF1'P(0,T; (L9(M3))3) p LP(0, T; {W^^'i{W^)f)] 

with the norm 



as well as 

with the norm 

We denote 
and 



|u|lwP'9(0,T) '■— l|u||H/i'P{0,r;L9(M3)) + ||u|| j;^p(o,T;W/2,9 

yj-.^ := n b];''^' X (tyi'^(M3))^° 



||(/,<7)||v'-:= 11/11 2(1-1) + 11/11 



n;(o, T) = ^^'"(0, T) n w2'2(o, r) 



Our first result is the following local existence and uniqueness: 

Theorem 3.1. Let Tq > be given and {uq, pQ, Eq) G Vq with p £ [2,oo),g € (3,cxd). 
There exists a positive constant do < 1, depending on Tq, p,, and \, such that if 

||(uo,po-l,^o)||yo < -^o, (3.1) 

then the initial-value problem (|1.3p - (|1.4|) has a unique strong solution on "M? x (0,To), 
satisfying u € >V(0,To) and 

{p -1,E)£ {W^''f{0, To; L«(M^) n L\M.^)) n 1^(0, To; W^'^^iR^) n W^'^M.^)))^" . 

The solutions in Theorem 13.11 is local in time since 6o = So{To) implies that Tq is finite 
for a given 6o <^ 1- 
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Remark 3.1. Notice that if q > 3, then by Theorem 5.15 in [Ij, the imbedding W^''^{ 
C^(M^) is continuous. Here, the notation C^(M'^) means the spaces of bounded, continuous 
functions in M^. Hence the condition (j3.ip imphes that, if we choose Sq sufficiently smah, 
by Sobolev's imbedding theorem, there exists a positive constant Cq such that 

po > Co > 0, for a.e. x G M^. (3.2) 

Remark 3.2. An interesting case is the case q < p. Indeed, by the real interpolation 
method and Theorem 6.4.4 in [3], we have 

Then, if we replace the functional space Vq''^ in Theorem 13.11 by 

V^'" := ((T^^(^-p)'^(M3))3n(^i-f''?(M3))3^ X (Tyi'^(R3))io, 
Theorem 13.11 is still valid. 



For the solutions claimed in Theorem 13.11 the following estimates hold uniformly in 
time. 

Theorem 3.2. Let be the same as in Theorem \3.1\ and assume < 5 ^ min{^,(5o}. 

(I) // the initial data satisfies ||(uo, po ~ li -E^o)IIV() ^ '^^) 0'''^^ fj, and A satisfy the assump- 
tion (16.18P below, then the solution {p, u, E) constructed in Theorem \3.1\ satisfies 

\Np\\l^{0,To;Li{R3)) < Cy^, ||Vi?||ip(o,To;L9(K3)) < Cv%- 

(II) // in addition, we assume p = 2 and q G (3,6]; and the initial data satisfies the 
compatibility condition (j2.2ip and 

dw{poFj) = 0, (3.3) 

/ (l:Po\vLo\^ + lpo\Eo\^ + ^—ipl-lPo + l-l))dx<5\ (3.4) 
Jrs \^ 2 7-1 J 

||Vuo||l2 + ||uo • Vuo||l2 + ||Auo||l2 + ||Vpo||l2 + llV-EolUanL'? < (3-5) 
Then, the solution {p, u, E) has the following improved estimates: 

[maxjg[o,T] max {11/5 - l||^yi,2n^yi,,(t), \\V ph^nLiit), \\E\\wi,2nw^'iii)} < C5o, 

< 3-9 3-9 

i I|VpIIl2(o,To;L'?(R3)) < C6o, ||5tu||i2(o,To;L'J(R3)) < C6q'^ , || Vii^||2,2(o,To;L'J(R3)) < C6q^ , 

for some 9 G (^, 1], where the constant C is independent of the time and 6. 

Remark 3.3. Under assumption (j3.3p . the authors in [151 [H] showed that the property 
will insist in time, that is, for all t > 0, div(yoF^) = 0. 

Remark 3.4. It is remarkable to point out that if we could obtain one better estimate on 
ll^/'llL2(o,r;L9(R3))) say p\\l'^(o^T;Li{r'^)) ^ C^o"*"" for some a > 0, then one could extend 
the local existence in Theorem 13.11 to the global existence for t £ [0,oo). 
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4. Local Existence 

In this section, we prove the local existence of strong solution in Theorem 13.11 using a 
fixed-point argument. To this end, we introduce the following new variables by scaling 

s := uh, y := ux, v{y, s) := ^u(x, t), r{y, s) := p{x, t), G{y, s) := E{x, t), 

where > will be determined later. Then, system (|1.3|) becomes 

rt + div(r't;) = 0, (4.1a) 

{rv)t + div {rv v) - fiAv - {fi + A)Vdivz; + v^'^VP = z^^^div [r{I + G){I + G)^^ , 

(4.1b) 

Gt + v-VG = VvG + Vv,. (4.1c) 
This scaling is needed in order to apply a fixed-point argument. From (|2.20p . one has 

VfcGij + GikViGij = VjGik + GnjVnGik- (4.2) 

Thus, if we denote by Gi the i-th row of the matrix G (or the i-th component of the vector 
G), then becomes 

curl Gi = GnjVnGik - GikViGij. (4.3) 

The proof of local existence of strong solutions with small initial data will be carried out 
through three steps by using a fixed point theorem. Instead of working on (|1.3p directly, 
we will work on (j4.ip . We note that (|4.ip is just a scaling version of (jl.3p . It can be seen 
from the argument below that we only need to verify the local existence in W^''^{0,T), 
< T < Tq, while the initial data belongs to Vq'"^- 

4.1. Solvability of the density with a fixed velocity. Let Aj{x,t), j = l,...,n, be 

symmetric mxm matrices in M" x (0,T), B{x,t) an m x m matrix, f{x,t) and Vo(x) be 
m-dimensional vector functions defined in M" x (0,T) and M", respectively. 
For the following initial-value problem: 

n 

dtV + M^^ t)djV + B{x, t)V = fix, t), 

j=i (4-4) 

V{x,0) = Vo{x), 
we have 

Lemma 4.1. Assume that 

A, E [C{0, T- H^{W^)) n Ci(0, T- H-'^-\W% "^"^ , j = 1, n, 

s e C7((o,r),F"-i(M"))™^", / G c((o,r),/7"(R"))'", t/q e ^"(k")"", 

with s > ^ + \ is an integer. Then there exists a unique solution to (14. 4p . i.e, a function 

V G [C([0,T),i7"(M"))nC7i((0,T),i7"-i(R"))]'" 
satisfying (j4.4p pointwise. 

Proof. This lemma is a direct consequence of Theorem 2.16 in [23] with Ao(x,t) = /. □ 
To solve the density with respect to the fixed velocity v E yV{0,T), we have 
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Lemma 4.2. Under the same conditions as Theorem \3.1{ there is a unique strictly positive 
function 

r := S{v) e W^'P{0, T; L^R^)) n L°°(0, T; ^^^'^(M^)) 

which satisfies the continuity equation ()4.1ap andr — 1 £ L°°(0, T; Moreover, the 
density satisfies 

l|Vr||icx,(o,T;L9(iR3)) < C{T, ||^^||w(o,T)) (l|Vro||ig(K3) + l) , (4.5) 
and the norm \\S{v) — l||p^/i,9(]g3)(t) is a continuous function in time. 

Here, and in what follows, C stands for a generic positive constant, and in some case, 
we will specify its dependence on parameters by the notation C(-). 

Proof. For the proof of the first part of this lemma, we refer the reader to Theorem 9.3 in 
|24j . or the first part of the proof for Lemma 14.31 below. The positivity of density follows 



directly from the observations: by writing (|4.1ap along characteristics ^X{t) = v, 

^r{t, X{t)) = -r{t, X{t))diYv{t, X{t)), X{0) = x, 
at 

and with the help of Gronwall's inequality, 

(infpo)exp j \\divv{t)\\Lao^^3)dx^ < r{t,x) < (sup/?o)exp \\dwv{t)\\Loo^^3)dx 

Now, we can assume that the continuity equation holds pointwise in the following form: 

dtr + rdivTj + v ■ Vr = 0. 

Taking the gradient in both sides of the above identity, multiplying by |Vr|'^~^Vr and 
then integrating over R^, we get, by Young's inequality 



1 d_ 

q dt 



Vr||^^.„3. < / |Vr|^|divt;|(i2; + / r|Vr|«"^| Vdivz;|(ia 



+ I \Vv\\Vr\''dx -- I vVlVrl'^dx 



^ ll^'^llLcRS) (I|Vv||l°o(r3) + ||r||ioo(K3)||Vdivt;||i,(R3)) (^-6) 

+ - I div^;|Vr|''dx + ||r||icx.(iR3)||Vdivf ||j;^9(K3) 
< C\\Vr\\lg,„3A\v\\w2,g/^3) + ||r||ioom3)||Vdiv't;||i,(jj3), 



since q > 3. Then (14. Sp follows from Gronwall's inequality. 

Finally, noting from ([i^ and (ji3]) that ^||Vr||^,.jg3. G L^(0,T), and hence 



|l|V(r-l)||L(^3)eLi(0,T), 

which together with ()4.5p implies that ||V(r — l)||^,^jg3j(i) is continuous in time, and hence, 
||V(r — l)||j;^ii(lg3)(i) is continuous in time. Similarly, from the continuity equation, we know 
that 

dt{r - 1) = -div((r - l)v) - divv G LP{0,T; L'^{M.^)), 
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which, together with the fact r - 1 E L'^{0,T; L'^iM.^)), yields r - 1 G C([0, T]; L9(M3)). 
Hence, the quantity ||r — l||^i,q(]K3)(t) is continuous in time. The proof of Lemma [^2] is 
complete. □ 

4.2. Solvability of the deformation gradient with a fixed velocity. Due to the 

hyperbolic structure of (j4.1cp , we can apply Lemma 14.11 again to solve the deformation 
gradient G in terms of the given velocity v G W(0,T). For this purpose, we have 

Lemma 4.3. Under the same conditions as Theorem \3.1\ there is a unique function 

G ■= T{v) e Tyi'P(0, T; L'^iM.^)) D L~(0, T; W^^'^iM.^)) 
which satisfies the equation ()4.1cp . Moreover, the deformation gradient satisfies 

l|VG||ioo(o,T;L.{M^)) < C(T, ||7;||w(o,T)) (l|VG(0)||i,(K3) + l) , (4.7) 
and, the norm \\G\\i^i,q(^3^{t) is a continuous function in time. 

Proof First, we assume that v G C^(0, T; ^^(M^)), Go G G^(M^). Then, we can rewrite 
(|4.1cp in the component form as 

dtGk + V ■ VGfc = V^; Gk + Vvk, l<k<3. 

Applying Lemma HT] with Aj(x,t) = Vj{x,t)I for 1 < j < 3, B{x,t) = Vv, and f{x,t) = 
Vvk, we get a solution 

oo 

G G Pi {G^(0,r,F'-i(M3)) n (0,r;/7'(M3))| , 
1=1 

which implies, by the Sobolev imbedding theorem, 

oo 

Gg f]C\0,T;C''{M.^)) = G\0,T;G'^{R^)). (4.8) 

k=l 

Next, for v G W{0,T), there are two sequences: Vn G G'^{0,T;C^{R^)), Gf^ G Go°°(M^), 
such that Vn^vin W{0, T), G[J ^ Go in ^^'^(R^), thus Vn ^ v in G{B{0, a) x (0, T)) for 
all a > where B{0, a) denotes the ball with radius a and centered at the origin. According 
to the previous result (j4.8p . there are a sequence of functions {Gn}'^=i C G^(0, T; G°°(M^)) 
satisfying 

dtGn + Vn-VGn = VVnGn + VVn, (4.9) 

with G„(0) = Gg. Multiplying (gj]) by \Gn\''~^Gn, and integrating over R^, using inte- 
gration by parts and Young's inequality, we obtain, 

1 d 
q dt 



\Gn\''dx 

= -- I Vn-y\GnVdx+ [ VVn\Gnr^Gldx+ [ VVnlGn^^Gndx 
Q 7R3 7ir3 J^3 

< — ^||G„||^<,(||Vi;„||loo + ||Vi;„||l9) + llVunlU?. 
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From Gronwall's inequality, one obtains, 

\Gn\Hx 



< 



< 

Thus, 



J ^ \Gn{OWdx + q llVvnllL.exp ^- ^ (g + 1)(||Vw„||l^ + \\Vvn\\Li)dT^ dsj 
X exp {q + l)(||Vt'„||Loo + \\Vvn\\Li)ds^ 

/ \Gn(,0)\''dx + q \\Vvn\\Lids] ex.p [ / {q + l){\\VVn\\L°° + \\'^Vn\\Li)ds 

Jr» Jo J \Jo 



l|G'n||L°°(0,r;L'/(R3)) < C(T, \\v\\lp(0,T;W2,i{R3))) (11^(0)112,9(183) + l) < OO. (4.10) 

Hence, up to a subsequence, we can assume that the sequence {vn} was chosen so that 
Gn^G weak-* in L~(0, T; L''(R3)). 

Taking the gradient in both sides of ()4.9p . multiplying by |VG„|'^~^VG„ and then 
integrating over R^, we get, with the help of Holder's inequality and Young's inequality, 

-—Wwr \\i 

< / \VGn\'^\Vvn\dx+ [ |G„||VG„|''-^|VVi'„|(i3; 

+ 1 \VVn\\VGn\'^dx-- I VnV\VGnVdx+ I \VV Vn\\V Gn^'^ dx 

< / \VGn\'^\VVn\dx+ [ \Gn\\VGn\'^'^\VVVn\dx 

+ / \VVn\\VGn\'^dx+- [ \VVn\\VGn\'^dx+ [ \VVVn\\VGn\'^~^dx 

Jr3 q Jrs Jr3 

< G\\VGn\\''rn(Tf.3\\\Vn\\w'^,'i(R3) + dlGnll' .5,3,^ + 1) ||pi/2,5m3) 



(4.11) 



< C||VG„||^ .j,3Jt'„||;y2,9(iR3) + (||Gn||Lm3N + l)lbn|lw2>'J(R3), 



|^q(IR3)ll'yn||vK2''J(R3) + (ll^'nll^q 
since q > 3. Using Gronwall's inequality and (j4.10p . we conclude from (j4.1ip that 

||VG„||ioo(o,T;L'J(R3)) < G{T, \\Vn\\w{0,T)) ( || VG„ (0) || ^9 {R3) + l) , 

and hence, 

l|VG'||ioo(o,T;L9(R3)) < liminf ||VG„||ioo(o,T;L 



(4.12) 

<C(T,||z;||w(o,T))(l|VG(0)||i,(i,3) + l). 

Therefore, 

\\Gn\\L^{0,T;W^'i{R3)) < C{T, ||'u||lp{o,T;W'2,9(ir3)) , || G(0) || (R3)) < OO. 

Furthermore, since g > 3, we deduce G G L°°{Qx) and 

\\G\\l°°{Qt) ^ C"!^' II^IIlp(0,T;W'2>'J(r3)), ||G(0)||vj/i,9(]j3)) < OO. 
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Passing to the limit as n — )• oo in (j4.9p . we show that (j4.1cp holds at least in the sense of 
distributions. Therefore, dtG £ LP{0,T; L'}{R^)), then G £ W^^p{0,T; Li{R^)), and hence 
G G C([0,r];L'?(R3)). 

Finally, to show that the quantity ||G||yi/i,<7(]R3)(t) is continuous in time, it suffices to 
show that ||VG||2,9(]r3') is continuous in time. Indeed, from (j4.1ip . we know that 

d 

which, with (j4.12p . implies that ||VG||ig(R3) G C([0, T]). The proof of Lemma 14.31 is 
complete. □ 

4.3. Local existence via the fixed-point theorem. In order to solve locally system 
(j4.ip . we need to use the following fixed point theorem (cf. 1.4.11.6 in 



Theorem 4.1 (Tikhonov Theorem). Let M be a nonempty bounded closed convex subset 
of a separable reflexive Banach space X and let F : M ^ M be a weakly continuous 
mapping (i.e. if Xn G M,Xn x weakly in X, then F{xn) — t- F{x) weakly in X as well). 
Then F has at least one fixed point in M . 

Now, let us consider the following operator 

Luj := ujt - ^Aw - (/i + A)Vdiva;, co E W^'^iO, T). 

One has the following theorem by the maximal regularity of parabolic equations; see 
Theorem 9.2 in [24j, or equivalently Theorem 4.10.7 and Remark 4.10.9 in [2] (page 188). 

Theorem 4.2. Given I < p < oo, ujq £ V^''^ and f G LP{<d,T;Li{m?f ), the Cauchy 
problem 

Luj = f, tG(0,r); uj{0)=uJo, 
has a unique solution uj := L^^iojQ, f) G W^'''(0, T), and 

ll^ll>Vp.9(o,T) < C (^||/||lp(o,T;L9(r3)) + ll'^o || yp--?) , 

where G is independent of ujq, f and T. Moreover, there exists a positive constant cq 
independent of f and T such that 



\\(^\\wp-i(o,T) > Co sup ||a;(t)||vP,9. 

te(o,T) " 

Notice that Theorem 14.21 implies that the operator L is invertible. Thus we define the 
operator n{v) : WP'i{0,T) ^ WP''i{0,T) by 

n{v) := L-^ (vo, dt{{l - S{v))v) - div(5(t;)t; C^v) + u~^V{P{l) - P{S{v))) 

+ i^-^diyiS{v)iI + Tiv)){I + Tiv)^)') . (4.13) 

Then, solving system (j4.ip is equivalent to solving 

v = 'Hiv). (4.14) 

To solve ()4.14p . we define 

BniO) := {v G WP'''iO,T) : \\v\\wp,,^o,T) < R}- 
Then, we prove first the following claim: 
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Lemma 4.4. There are z^, T > 0, and < R < 1 such that 

n{BR{o)) c Br{o). 

Proof. Let T>0, 0<i?<l and v G Bji{0). Since S{v) solves ()4.1ap . we can rewrite 
operator 7i as 

n{v) = L-^(vo, (1 - Siv))dtv - S{v)v • Vt; + u-'^V{P{l) - P{S{v))) 

+ i^-^dw{S{v)iI + Tiv)){I + T{v)y)) . (4.15) 



Thus, it suffices to prove that the terms in the above expression are smah in the norm of 
LP{0,T; (L''(M3))3). 

First of ah, we begin to deal with the first term by letting r := S{v) — 1. Thus, r 
satisfies the equations 

dtf + div(ru) = 0, r(x, 0) = ro — 1. 
Repeating the argument in Section 14.21 again, we obtain 

II^IIl°°(Qt) - C'll^llL°°(0,T;iyl.9(M3)) < C\\ro - 1 1| ,y l,9(]j3) C(r, ||f || >v(o,T) ) 

< ||ro-i||^i,,(K3)C(r,i?) <c(r)i?, 

where, by the formula of change of variables, we deduce that 

3 3^ 

Iko - 1||l9(r3) < ui \\po - 1||l9(]r3) < R, ||Vro||L,(K3) < u i ||Vpo||m(r3) < R, 



if \\po — 1||l9(ir3) is small enough and > 1 is large enough. Hence, due to the assumption 
v S Br{0), we obtain 

II (1 - S{v))dtvhr>io,T;L'>(m) ^ C{T)R\ (4.16) 

Secondly, by the Sobolev imbedding, 

\vVv\''dx < \\v\\\q,^3^\^v\\\^,^z. < C||u||^,.]g3J|f ||^2,gm3), 



and thus, since Vri'P(0, T; L'?(M3)) ^ C([0, T]; ^^(M^)), we deduce 
Therefore, we get 

\\S{V){V ■ V)v\\LPio,T;LHm) ^ (^-l^) 

Thirdly, for the term VP{S{v)), we can estimate it as follows 
||VP(5(t;))|| 

< C(r) sup {P'(r/) :C(T)-i <r/<C(r)}(||Vro|U,(M3) + l). '^^^ 
Fourthly, for the term div(5(t;)(/ + T{v)){I + T(v))^ , we have 
|div(5(z;)(/ + Tiv)){I + T{v)y)\ < |V5(7;)||/ + Tiv)\^ + 2S{v)\VT{v)\\I + Tiv)\, 
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and hence, 

\\div{S{v){I + T{v)){I + r(t;))T)||iP(o,T;L.) 

< \\\VS{v)\\I + T{v)\^\\Lvio,T-M) + nS{v)\VT{v)\\I + T{v)\\\lv(o,t-m) (4-19) 

< C{T)M, 

with 

M = max{||Go||vKi>'j + I|Go||l°°(ir3) + 1^ Ikollwi^-? + 1, lko||L°°(R3) + l}^ < oo. 

Combining together (I4.16p -( j4.19p . using the Theorem 14.21 and assuming parameter v 
sufficiently large and R <1 sufficiently small, we get 

\mv)\\y,i^o^T)<C{T){R^ + u-^)<R. 

The proof of Lemma 14.41 is complete. □ 

Thus, it is only left to show the following: 

Lemma 4.5. The operator % is weakly continuous from W'''(0, T) into itself. 

Proof. Assume that Vn ^ v weakly in >V^''^(0,T), and set r„ := S{vn), Gn ■= T{vn), then 
{rn}n=i and {Gn}^=i are uniformly bounded in L°^{0,T;W^''i{R^)) D ^^^'^(O, T; L9(R3)) 
by Lemmas l4.1l and l4.3i Hence, up to a subsequence, we can assume that rn ^ r and G„ — t- 
G weakly* in L'^{0,T-W^'i{R^)) n W^'P{0,T; Li{R^)) and then strongly in C((0,r) x 
5(0, a)) for all a > 0. And at least the same convergence holds for f„. Thus, (j4.1ap and 
(|4.1cp follow easily from above convergence. 

Since r„ ^ r weakly* in L°°{0,T; W^'''{M.^)) n W^'P{0,T; Li{M.^)), we can assume that 
P'(5(u„))V5(?;„) ^ P'{S{v))VS{v) weakly in LP{0,T; Li{R^)) and hence, 

(0, VP{S{vn))) L-^ (0, VP{S{v))) weakly in W(0, T), 

since the strong continuity of L^^ from L^(0, T; L''(M^)) into W(0, T) and the linearity of 
the operator L imply also the weak continuity in these spaces. 

Similarly, since dtVn dtv weakly in LP(0, T; L<?(M3)) and r„ r in C((0, T) x B{0, a)) 
for all a > 0, we have (rg — rn)dtVn — )• {re — r)dtv weakly in L*'(0, T; L''(]R^)) and conse- 
quently 

(0, (re - rn)dtVn) ^ (0, (r^ - r)9t^;) weakly in W(0, T). 

Since V?;„ Vv weakly in W^1'P(0, T; Ty-i'''(M3)) n LP(0, T; ^"^'''(M^)) which is com- 
pactly imbedded in to C{[0,T]; L''(B{0,a))) for all a > 0, we can assume that Vn ^ v 
strongly in L°°(0, T; L9(S(0, a))) for ah a > 0, and then 

weakly in LP(0, T; L'?(m3)). Hence 

L-^ (0, S{vn)ivn ■ V)vn) ^ L'^ (0, 5(v)(t; • V)v) weakly in >V(0, T). 

Finally, due to the facts that r„ ^ r and G„ ^ G weakly* in L°°{0,T;W^''^{R^)) D 
W^'P{0,T;L''{M.^)) and strongly in C((0,T) x 5(0, a)) for ah a > 0, we deduce that 

dw{Sivn){I + T{vn)){I + Tivn)V) ^ div(5(^;)(/ + T{v)){I + T{v))^) 
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weakly in LP(0, T; L'?(M3)). Therefore, 

(O, dw{S{Vn){I + T{Vn)){I + r{Vn)V)) ^ L'^ (o, dw{S{v){I + T{v)){I + T{v))^) 



weakly in >V(0,T). 

Thus, we can conclude that 

n{vn)^n{v) weakly in >V(0,r), 

due to the weak continuity of map L^^(v,0). The proof of Lemma 14.51 is complete. □ 

Therefore, by Theorem 14. H there exists at least a fixed point 

V = n{v) G Br{Q) C W(0, T), (4.20) 

of mapping %. The fixed point v provides a local in time solution (p, u, E) of system (jl.3p 
near its equilibrium through the scaling with v sufficiently large. 

The proof of the local existence in Theorem 13.11 is complete. The uniqueness will be 
proved in the next section. 

5. Uniqueness 

In this section, we prove the uniqueness of the local solution found in the previous 
section. Notice that, the argument in Section H] yields that 

e L2(0,r;L2(R3), Vr e L2(0,r;L2(M3)), VG e L^{Q,T-L^{W^)). 

Hence, using the interpolation, we deduce that 

G LPo(0,r;L3(M3), Vr G LP"(0,r;L3(M3)), VG e L'p°{{),T-L^{W^)), 

where 



po 2 p 3 2 q 
for some 9 G [0, 1]. Now, assume that vi, V2 satisfy (|4.20p for some T > 0. Let 

r := S{vi) - S{v2), v:= vi - V2, G := T{vi) - T{v2), 

with a little abuse of notations (however, there should be no confusion in the rest of this 
section). Then, we have 

dtr + vi ■ Vr + v ■ VS{v2) + rdivfi + 5(f2)divf = 0, r(0) = 0. (5.1) 

Multiplying (j5.ip by r, and integrating over M^, we get 

1 d 



2dt 
which yields 
d 



- [ \r\'^divvidx + [ ( vV S {v2)r + \r\'^ divvi + rS {v2)divv) dx = 0, 
2 Jr3 7k3 



^^ii'^lli2(K3) ^ l|divui||Loc||r||^2 +e||VT;||^2 + C(e)||V5(t>2)r||^6 



+ e\\Vv\\i,^^s) + CXe)\\Siv2)\\io.\\r\\i 
< \\divvi\\Loo\\r\\l2 + e\\Vv\\l2 + C{e)\\VSiv2)\\l3\\r\\l2 (5-2) 

+ e\\Vv\\l,(^^,^ + C{e)\\Siv2)\\loo\\rg 
<m{e)\\r\\l2+2e\\Vv\\l2 
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where e > 0, ??i(e) = ||div7;i||Loo + C{e) {\\V S {v2)\\l, + ||5(t;2)||ioo) . 
Similarly, from ()4.1cp . we obtain 

dtG + vi-VG + vVG2 = VviG + VvG2 + ^v, G(0) = 0. (5.3) 

Multiplying (j5.3p by G, and integrating over M^, we get 

1 d 



2dr 



G\\h --I \G\^d\wvidx + [ V- VT{v2) : Gdx 

\G\'^Vvidx + / VvT{v2) : Gdx + [ Vv: Gdx, 
;3 Jr3 J^3 



which yields 

j^\\G\\l.(^^,)<\\divvi\\L^\\G\\l,^ 

+ 6||Vt-||i.(j,3) + C(£)||r(^;2)||ioo||G||2 + e||V^||i.+C(e)||G||i. 
<||divr;i||^^||G||i.+e||Vt;||i.+C(e)||Vr(t;2)|li3||G||i. 

+ ^l|V^|li2(K3) + C(e)||r(^;2)||io.||G||i + e||V^||i.+C(£)||G||i. 

<%(e)l|G||i2 +3e||Vt;||i2(K3), 

where 772(e) = ||divt;i||ioo + C{e) (||Vr(7;2)|li3 + ll'r(7;2)|lioc + l) . 
For each vj, j = 1,2, we deduce from (|4.1b[) that 

S{vj)dtVj — fiAvj — + X)Vdivvj 

= -S{vj){vj ■ V)vj - VP{S{vj)) + div(5(7;,-)(^ + T{vjyKl + Tivj)^), 
Vj{0) = vo. 

Subtracting these equations, we obtain, 

S{vi)dtVi — S{v2)dtV2 — fJ.Av — (^ + A)Vdivz; 

= -S{vi){vi ■ V)t;i + S{V2){V2 ■ V)V2 - VP{S{vi)) + VP{S{v2)) (5.5) 
+ div(cS(^;i)(/ + T{vi)){I + T{vi))'^) - dw{S{v2){I + T{V2)){I + T{v2)y). 

Since 

— S{vi){vi ■ V)vi +S{V2){V2 ■ V)V2 

= —S{vi){v ■ V)vi — iS{vi) — S{V2)){V2 ■ \^)vi — S{v2)iv2 ■ V)^;, 

and 

S{vi){I + T{vi)){I + T{vi))'^ - S{v2){I + T{v2)){I + T{v2)V 
= S{vi)G{I + T{vi))~^ + r{I + T{v2)){I + r(fi))"^ + S{v2){I + T{v2))G^, 
we can rewrite (|5.5p as 
<S{vi)dtv — jsAv — (/u + A)Vdivt; 

= -rdtV2 - S{vi){v ■ V)vi - S{v2 ■ V)vi - S{v2){v2 ■ V)v - VP{S{vi)) + VP{S{v2)) 
+ dwiS{vi)G{I + T{vi))^ + r{I + T{V2)){I + nvi)Y + S{V2){I + Tiv2))G'^). 

(5.6) 
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Multiplying (j5.6p by v, using the continuity equation (jl.3ap and integrating over M^, we 
deduce that 

~l S{vi)\v\^dx+ [ UVv\^ + {fi + X)\divv\^)dx 
^ at Jm3 Jr3 

-S{vi){vi ■ V)v ■ V — rdtV2V — S{vi){v ■ V)viv — S{v2 ■ V)viv 



- S{v2)iv2 ■ V)vv - VP{S{vi))v + VP{S{V2))V 

- {S{vi)G{I + Tivi))'^ + r(I + T{V2W + T{vi))^ + S{V2){I + T{v2))G'^)Vvdx 
e\\Vv\\l,+C{e)\\S{v^)\\l^\\v,\\l^^^^^^^ 

+ ||5(?;i)||Loo||V?;i||L^||t;||22 +2||z;2||L-||V7;i||Loo(||r||22 + ll^;!!^^) 
+ £||VHli2+C(e)||5(t;2)||ioo||^2||ioo||Hli2+e||Vz;||i. 
+ C(e)(sup{P'(7?) : C(r)-i < r? < C{T)}f\\r\\l, + e||V^||i. 
+ C(e)(||5(^i)||i.(l + ||r(t;i)||ioo)||G||i. 

+ ||5(^;2)||ioo(l + ||r(^;2)||ioo)||G||i. + ||r||i.(l + ||r(t;i)||ioo)(l + \\T(v2)\\U)) 

2 



(5.7) 



<5e\\Vv\\i,+r^,ie)i\\r\\i, + \\v\\i, + \\G\\i^ 
with 

= C(e)||5K)||i^|bi||i^ +C(e)||att;2|li3 + \\S{v^)\M\Vvi\\l^ 
+ 2||?;2||l-||V?;i||l^ + C{e)\\S{v2)\\loo\\v2\\l^ 
+ C(e)(sup{P'(7?) : C{T)-' < r? < ^(r)})^ 

+ C{e)mv,)\\l^{l + \\nv,)\\U) + + ||r(^2)||ioo) 

+ (i + lin^i)llioo)(i + l|rMllioo)). 

Summing up ()5.2p . (|5.4p . and (j5.7p . by taking e = we obtain 

^ / (5(i;i)|?;|2 + |r|2 + |G|2)dx+ / (f,\Vv\^ + (fi + X)\divv\'') dx 
dt Jm3 Jk3 

< 2(r?3(e) + r?2(e) + + Mh + \\G\\h) (5-8) 

<2r]{£,t) [ {S{vi)\v\^ + \r\^ + \G\^)dx, 



with 



min{min3,gR3 5(t>i)(x, t), 1} ' 

It is a routine matter to establish the integrability with respect to t of the function r/(e, t) 
on the interval (0,T). This is a consequence of the regularity of vi,V2 € W(0,T) and the 
estimates in Lemmas 14. 21 and 14. 31 for S{vi), T{vi) with i = 1,2. Therefore, ()5.8p . combining 
with Gronwall's inequality, implies 



/ 



(5(fi)|t;|2 + |r|2 + = 0, for ah tG(0,T), (5.9) 
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and consequently v = 0, r = 0, G = 0. Thus, the uniqueness m Theorem 13. H is estabhshed. 



6. Uniform A Priori Estimates 



Up to now, we prove that for any given Tq, we can find a unique solution to the scaling 
system ()4.ip . That is, we have proved the local existence of solution to the viscoelastic fluid 
system ()1.3p and its uniqueness. For the unique solution we constructed in the previous 
sections, we have some a priori estimates uniform in time as stated in Theorem I3.2[ In 
this section, we prove the first part of Theorem 13.21 To simplify the presentation, we will 
focus on the case v = I, that is, system (jl.3p . 

We introduce the new variable: 

(T := V In p. 

Then, we have 



Lemma 6.1. Function a satisfies 

dta + V(u • a) = 0, 

in the sense of distributions. Moreover, the norm ||o"(t)||L9 



(6.1) 



is continuous in time. 

Proof. We follow the argument in |24) (Section 9.8) by denoting as = Ssa, where is the 
standard mollifier in the spatial variables. Then, we have 



with 

Ue = V(U • ae) - 

=: 7^i +7^?. 



dtae + V(u • as) = Ue, 
SeV{u ■a) = {u- Vae - ^^(u • Va)) + {a^Vu - Se{a • Vu)) 



3.2) 



Since a G L°°{0,T; L'i{R^)) and u G LP{0,T;W^^°°{R^)), we deduce from Lemma 6.7 in 
(cf. Lemma 2.3 in [T9]) that U] ^ as e — )• 0. Moreover, 



{as - o-)Vu||ii(o,7..i, 



< Ik 



iVul 



Lp(0,T;L°°(R3)) 



and Se{a • Vu) ^ cj • Vu in L\0,T; Li{R^)) since a • Vu G LP{0,T; L'^{R^)). Thus, we 
have 7^2 ^ in LP{0, T; L'}{R^)). Then, taking the limit as e in ([621), we get ([eT]) . 
Multiplying ()6.ip by |cr|''~^cr, and integrating over R^, we get 



d 



dt 



LI 



f {—djUj.ajak\a\'^ ^ divu|(T|'')(ix 

JR3 q 



1 



< ||Vu||Loo||fT||^, + -||divu||ioo ||cr|||, < C||u|||y2,9||cr|||^g. 



Dividing the above inequality by ||cr||^, ^, we obtain 

d „ , 



dt 



\0'\\Li 



< C\\u\\iY2,q \\a\\Li. 



Since a G L°°{0,T;Li 
Lemma l6. II is complete. 



, jt\\a\\L^ G LP{0,T). Thus, \\a\\L^ G C{0,T). The proof of 

□ 
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For a given R = 60 <^ 1 as in Section^ if the initial data satisfies ||u(0), po — 1, -E'(0)||vb — 
6^ with < (5 <C min{^, Jo}, let T{R) be the maximal time T such that there is a solution 
of the equation u = 'H(u) in i?/j(0). By virtue of Lemma [4.2[ Lemma [4.31 and Lemma EH 
we know that ||5(u) — l||yi/i,q(]K3), ||<t||l9 and ||T(u)||54/i,q are continuous in the interval 
[0,T(R)). On the other hand, under the assumptions on initial data and Remark 13.11 we 
know, if 6 is sufficiently small, then 

||cT(0)||i,(K3) < 7^||Vp(0)|U,(K3) < 55 < L 

Hence, there exists a maximum positive number Ti such that 

max{\\S{u) -l\\wi,,{t),\\a\\L^{t),\\T{u)\\wi,,{t)} <VR<^1 for ah tG[0,ri]. (6.3) 

Now, we denote T = min{T(i?), Ti}. Without loss of generality, we assume that T < 00. 
Since g > 3, we have 

\\p - i|Il°°(r3) < c\\p - i||vi/i,9(]K3) < cVr < - 



2' 



if R is sufficiently small. Hence, one obtains 



1 3 
- < p< -. 



On the other hand, for any given t £ (0, T), we can write 
l|u(t)||i, = ||u(0)||^,+^*^||u(.)||^,d. 

= iiuoiii. + - r (Mm^' [ \u{s)rMs)dsu{t)dx) dt 

Q Jo \ ilR3 / 



Q Jo 



<I|uo|Il9 + - / l|u(s)||i/||5sU||Lg(is 



t ^—^ t - 



< s'^P + ?-RP, 

q 



and consequently. 



l|u||Loc(o,t;i,) < i^S^P + ^RPj ' <CR, te (0, T). (6.4) 
Similarly, we have, for all t £ [0,T], ||u||2,oo(o,t;L2) < CR. 

6.1. Dissipation of the deformation gradient. The main difficulty of the proof of 
Theorem 13.21 is to obtain estimates on the dissipation of the deformation gradient and the 
gradient of the density. This is partly because of the transport structure of equation (ll.3cp . 
It is worthy of pointing out that it is extremely difficult to directly deduce the dissipation 
of the deformation gradient. Fortunately, for the viscoelastic fluids system ()1.3p . as we 
can see in [U [I3l [Ml [TSj [I6l [iTl [18] , some sort of combinations between the gradient of 
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the velocity and the deformation gradient indeed induce good dissipation. To make this 
statement more precise, we rewritten the momentum equation ()1.3bp as, using (|1.3ap 

dtu - ^Au + A)Vdivu - divE = -p{u ■ V)u - VP{p) + div(/)(/ + E^) 

+ dw{{p - 1)E) + divipEE'^) + (1 - p)dtu, ^^'^^ 

and prove the following estimate: 
Lemma 6.2. 

I|V£^||lp(o,T;L9(r3)) < C{p, q, p) (^R + VR\\cr\\LP{o,T;Li(^^))) ' (6-6) 
Proof. Now we introduce an operator £ : w i— Coj such that 

-pAC{Lo) - {p + A)Vdiv£(cj) = u. 

Then we denote Zi{x,t) as 

Zi:=C{dwE). (6.7) 

Notice that, from the elhptic theory, E £ LP{R^) with 1 < _p < oo, then Zi £ W^'P 
Then, (j6.5p becomes 



dtu - /xA (u - Zi) - (/i + A)Vdiv(u - Zi) = (6.8) 
where, with the help of Remark 13. 3^ 

J^i = -p{u ■ V)u - VP{p) + div{{p - 1)E) + divipEE'^) + (1 - p)dtu. 
Also, from (jl.3p . we have 

_ / d(^\ ^^(^^(dE\\^ ^^^^^^ ^ _ _ ^^^^^ _ ^^^g^ 



at \ dt J \ \dt 
From (|6.8p and (|6.9p . we deduce, denoting Z = u — Zi, 

dtZ - pAZ -{p + A)VdivZ = F ■.= Fi- F2, (6.10) 
where T2 = C (div(Vu + Vu£' — (u • V)E')) . Equation (I6.10p with Theorem 14.21 implies 

||Z||w(0,T) < C{p,q,p,X) ^||Z(0)||^2(i_i) + ||-^||LP(0,r;L9(R3))^ 

< C{p,q,p,X) [R+ ||-7^||lp(o,T;L'j(r3))) . 
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Next, we estimate ||-^j||Lp(o,T;L9(R3))) « = 1,2. Indeed, for J^i, using (j6.4p . we have 

I|-^i||lp(0,T;L'J(M3)) 

< l|p||L°°(QT)ll"lli-°°(0,T;L'J(R3))||Vu||iP(0,T;L°°(R3)) 

+ "II V£'||LP(o,T;L'3(f7)) + a|k||LP(o,T;L'3(f7)) II -^^11 L°° (Qt) 
+ \\P - 1||l°°(Qt)II^-^IUp(0,T;L'?(R3)) 

+ lkllLP(0,T;L'?(R3))||£'||Loo(Qy) + ||o-|| lp(o,T;L'J(K3)) |cx) (q^,) 

+ ll-^llL°°{QT)ll^^llLf(0,T;L'J(R3)) + ||/3 - 1 1| L°o (Q-p) || ||iP(o,T;L9{R3)) ^^-^^^ 

< 2i?||u||^p(o^T;W2>'J(]R3)) + a||V£'||/,p(o,T;L'3(R3)) 

+ "lkllLP(0,T;L9{R3))||-E||Loo(Q^) + \/^||o-||lp(0,T;L'J(R3)) 
+ -R|kllLP{0,T;L9(R3)) + ^/R\\^ E\\i^v(Q^T;Li{M?)) + -^^ 

< i^a + a||V£'||ip(o,T;L9(K3)) + ^^I|o"IIlp(0,T;L'J(k3)) + V^II V-EH iP(o,T;L9(R3)) . 

Here, a = sup {a;P'(x) : ^ < x < |} and in the first inequality, we used the identity 

Vp = -pdivE^ - VpE^ 
due to Remark 13.31 And, for F2, we have 
||Vu + Vui?-(u.V)i?|| ^ 

< 



|u||LP(0,r;W'2.9n^2) + ||Vu||ip(Q ^.i3)||£'||ioo(o^r;L9) + l|u||LP(o,r;L3)||V£'||ioo(o,r;L9) 



<R + m. 

Here, we used the following Gagliardo-Nirenberg inequality 

T ~ 
l|Vu|| J, < ( / (||u(s)f||Au(s)||^-^)''dsy < ||u||ip(o,T;W2,,njy2), 

with 6* = T^zg- Hence, one can estimate, by LP estimates of elliptic operators, 

ll-7^2||LP(0,T;LnM3)) < ||-^2 H^^^^^^^^i,^^ 

<C(/i,A)||Vu + Vui?-(u.V)i?|| 3, (6.13) 

< C(/x,A)(i? + i?i). 



Therefore, from (j6.12p and (|6.13p . we obtain 

II-^IIlp(0,T;L9(M3)) < C(/i, A)(i?2 + i?) + a||V£'||ip(o,T;L9) + ^^I|o"IIlp(0,T;L9) 
+ ^/i?||VS|Up(0,r;L.). 

Inequalities (|6.1ip and (|6.14p imply that 

II^IIlp(0,T;VK2,9(ir3)) 

< C(p, g, ^, A) (i? + a||V£;||LP(o,T;M) + V^lkllLP(0,T;L«) + ^^^11 V^||lp(o,T;L9) 



(6.14) 



(6.15) 
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Hence, we have, from (j6.7p 

||div£J||2,p(o,T;L9(R3)) < f^\\Z\\LP{0,T;W^'i(R»)) + W'^W LP {0,T;W^-<1{ 
< C{p, q, H, A) (^R + VR\\cr\\LP{0,T;Li{R3)) + ^^11 V£;||LP(o,r;L9(R3))^ (6.16) 

+ C{p, q, /i, A)/ia|| VS||ip(o,T;L'?(]R3)). 
On the other hand, from the identity ()4.3p . we deduce that 

||curl £^i||LP(o,T;L'J(R3)) < 2 ||i?|| ^oo (g^,) || Vii^|| 2,p(o,T;L9( 



^ (6.17) 
Combining together (j6.16p and (j6.17p . we obtain 

I|V£^IIlp(0,T;L'J{R3)) < C{p, q, fl, A) (^R + VR\\cr\\LP{0,T;Li) + V^ll V-E||ip(o,T;L'3)) 
+ (^(P, Q, fJ', A)Ma||V^||ip(o,T;L'/]), 

and hence, by choosing \/R ^ ^ and using the assumption 

C{p,q, p, X)fia < 1, (6.18) 
one obtains (|6.6p . The proof of Lemma 16.21 is complete. □ 

Remark 6.1. The assumption (j6.18p is reasonable, because if we consider the special 
case: < p < 1 and A = 0. Then, after the scaling, we get a control on the constant 

_ 6+2 2q-6 

C{p,q, p, X) < C{p,q)fi ^1 , and hence C{p,q, p, X)fia < ap 3? — )- 0, as /i — )• 0. 



Remark 6.2. Notice that, in view of the above argument, estimate (16. 6p is actually valid 
for aU t £ [0,r], that is, for ah t £ [0,T], 

ll^-^llLP(0,t;L9(R3)) < C{p,q,fJ.,X) (R + V^||cj||LP(o,t;L9(R3 



6.2. Dissipation of the gradient of the density. To make Theorem 13.21 valid . we need 
further the uniform estimate on the dissipation of the gradient of the density. 

Lemma 6.3. For any t G (0,T), 

lk||LP(o,t;L9(R3)) <C{p,q,p)R. (6.19) 

Proof. Multiplying (fObj) by (t|(t|'^ ^ and integrating over M^, we obtain 
p + X d 



yR3 

= / (//Au + (/x + A)Vdivu) • cj|a|''"^(ix - / pdtu ■ a\a\'^~'^ dx 
- / p{u-V)u-a\a\'^~^dx-{n + X) [ V {u ■ a) ■ a\a\'^'^ dx 
+ / div{p{I + E){I + E)'^)-a\a\'^-^dx. 

JR3 

We estimate the right-hand side of (|6.20p term by term. 



(6.20) 



/ (/xAu + (;u + A)Vdivu) • cr|cr|''~^dx 

JR3 



^ 5-1 
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/ 



pdtVL ■ ct|cj|'' '^dx 



< \\dtu\\Li\\(T\\%'^; 



pu 



■ Vu • cr|cr|«"^(ix 



< ||u • Vu||l9||ct|||^/ < ||u||L^r||u||^2,5||o-|||^/; 



/ V {u ■ a) ■ a\a\'^ ^dx = / dj\ik(Tk'^j\(T\'^ "^dx + / Ukdjdk{ln p)dj In p\a\'^ ^dx 

Jr^ Jr3 Jir3 

/ djUkcrk<7j\a\'^~'^dx + [ Ukdk\(T\'^\cr\'^~'^dx 

/ 9jUfc(Tfc(Tj|(T|''~^(ix + / Ukdk\(T\\a\'^^^ dx 

Jr3 J^3 

I djUkCTkcrj\a\'^^^dx -\ — / Ukdk\(7\'^dx 
Jrs Q Jrs 

/ djUkcrkcrj\a\^~'^dx / |cT|'^divudx 

Jr3 q Jr3 

< C\\Vu\\Loc\\a\\lg < C\\u\\iY2,q\\a\\lq, 
and, due to ()2.9p . we can rewrite 

di.ipiI + E)iI + EV)) =W£,+Mfl±M 

/ i OXj 



dxi 



dxi 



p{ej + Ej 



dxj' 



then one has 

div(p(/ + E){I + E^) ■ a\a\'^-'^dx 
< \\VE\\L4l + E\\L^h\\%' < 2\\VEh.^^3)\\cT\\l-' 



I \ dEi \n — 1 

p\ej + Ej)— — ai\aY dx 

OXj 



On the other hand, we have 



pP\p) = P'(l) + {p-l) C {P'iviP - 1) + 1) + iviP - 1) + ^)P"{viP - 1) + 1)) dv, 

Jo 

for p close to 1 and consequently 

\\pP'{p)-p'{i)U^<\\p- 

<CVRsup!^\fix)\ : ^<x<^^ <CVR, 



1 3 
iLoo sup<! |/(x)| : - <x < - 
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where f{x) = P'{x) + xP"(x). Thus, from (|6.20p . we obtam 

^^|ikiii.+F'(i)ikir„ 

+ C\\u\\w2,,\\a\\L'i + C||V^||i,(K3) + VR\\a\\Li 

< C'lkllL9^(ll^u||L9 + ||9fU||L9 + ||u||^2,q||u||L9 

+ ||u||^2,5||cj||l9 + ||VS||i9(K3) + VR\\a\\Li 
and hence, by assuming that R <^ 1, one obtams 

W\\% + Ip'{1)\\c7 

< C'lkllL9^(ll^u||i9 + + I|u||w2,9||u||l9 + ||u||ty2,9||0-||L9 + ||VS||l9) 



q dt" "^^ 2 ^ ^" (6.21) 



Multiplymg dOTTl by || fll^g'^, we obtain 

fi + X d p 1 ^ II IIP 

^^IHIl. + 2^(i)IHIl. 

< C'lkllL9^(ll^u||L9 + + l|u|liy2,9||u||L9 + ||u||^y2,9||cr||L9 + ||V£'||l'3). 

Integrating the above inequahty over the interval (0, t), one obtains, by using (j6.6p . 
+ lp'(l)j\\a\\l,ds 



p 



1 

t 



+ (l|u||L°°(0,t;L'?(M3)) + lkllL°°(0,t;L'?(R3)) + l) ( / ||u||^2,9C^f 
+ l|V-^llLP(0,t;L'J(M3)) 

^llLP(0,r;L9(R3)) 

p 



+ C{p,q,n)R\^J^ \\a\\lgdsj (l + lkl|L°°(0,t;L'?(R3)) + l|u||L°°(0,t;L'J 

and hence, by letting R be so small such that C{p,q, fi)VR < one obtains 
H + X 



P ^ Jo 



p-i 

(6.22) 



1^ VJo 



26 XIANPENG HU AND DEHUA WANG 

Plugging (|6.4p into (|6.22p . we obtain 



p „CT(t)||^, + -P'(l) / |k||^„d5 



p-l 



/o 



<^||a(0)||^, + C(p,(?,/i)i?^VllI,d^) ' (l + lkllL-(o,i;M 
Then, Young's inequality yields 

P "JO 
< t^±AslP + C{p, q, fi)RP{l + MLo.^o,t;Li)Y, 



(6.23) 



for all < t < T. 

Now, we let R be so small that 



Due to the fact that ||o'(0)||j;^<j(jr3) < 6^, we can assume that ||o"(t)||L9 < in some 

maximal interval (0, tmax 

) C (0,r). If W < T, then, ||a(tmax)||L« = ^VR and by 
i/R= < Cip,q,p)pR{l + VR) < ^VR, 

which is a contradiction. Hence, imax = T and 

\\cr\\Li < ^VR, for all t G [0,T]. (6.24) 
Thus, by ()6.23p . one obtains ()6.19p . The proof of Lemma 16.31 is complete. □ 

We remark that, from ()6.6p and (j6.19p . one has 

l|V-E||ip(o,T;L9(R3)) < C{p,q,ij)R. (6.25) 

7. Refined Uniform Estimates 

In this section, we prove the second part and thus complete the proof of Theorem [37 
Define 

^max := sup|t > : 3 u G W(0,T) with u = 'H(u), such that, ||u||w(o,T) < R^ 
ll'5(u) - l||L°°(o,T;iyi.9) < lkllL°°(o,T;L9) < and 
ll'7'(u)||ioo(o,T;iyi.9) < ^/^}, 
where ii was constructed in the previous section. 
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7.1. Uniform estimates in time. We now establish some estimates which are uniform 
in time T. First we prove the following energy estimates: 

Lemma 7.1. Under the same assumptions as Theorem \3.1{ we have 

I|Vu||l2(0,T;L2(m3)) < CR^, (7.1) 
I|u||l°°(0,T;L2{R3)) < CR^, (7.2) 
I|-^IIl°°{0,T;L2(R3)) < CR^, (7.3) 

\\P - 1|Il°°(o,T;L2(m3)) < CR^, (7.4) 
where C is a constant independent ofTG (0, Tmax)- 
Proof. First we recall that 

u G W^'^{0, T; l2(m3)) p ^^(0, T; VF^.z^j^a^^ 

and 

Multiplying equation ()1.3bp by u, and integrating over M'^, we obtain, using the conserva- 
tion of mass (jl.Sap . 

d 



dt 



[ (lp\u\^ + ^—{p'y + j-l)]dx+[ (^|Vu|2 + (/i + A)|divu|2)dx 
Jr3 \2 7-1 J J^3 



pFF~^ : Vudx. 



Here, the notation A : B means the dot product between two matrices. Thus, we have 

/ fip|u|2 + ^-(pT + 7-l)^ dx+ / [ {p\Vu\^ + {p + X)\dwu\^)dxds 
Jr» \^ 7-1 / Jo Jm.3 

= 1 (l-po\uo\^ + ^—-{pl + j-l)]dx- I I pFF~^ -.Vudxds. 
JR3 V2 7 - 1 J Jo J^3 

From the conservation of mass (|1.3a|) , one has 

[ fip|u|2 + ^-(pT-7p + 7- 1)^ dx+ / [ {p\Vu\^ + {p + X)\diYu\^)dxds 
Jr3 7-1 / Jo iR3 

= / (^PoluopH 3— -(p2; - 7/)o + 7 - 1) - / / pFF^ iVudxds. 

Jr3 \^ 7-1 / Jo Jr3 

(7.5) 

On the other hand, due to equations (jl.Scp and (jl.Sap . we have 

Ft (^'^1') = I'^l' + 4 = I'^l' + ^ (^"^ - " • 

= ||F|^ + 2pF:(VuF)-pu.V|F|^ (T.6) 

= ^IFP + 2pF : (VuF) + div(pu)|Fp - div(/)u|FP) 
at 

= 2pF : (VuF) - div(pu|Fp). 
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Integrating (|7.6p over R'^, we arrive at 

ld_ 

Since 



/ p\F\'^dx= / pF:{VuF)dx. (7.7) 



I I pF : {VuF)dxds = [ [ pFF^ iVudxds, 
Jo Jr^ Jo Jr^ 

we finally obtain, by summing (|7.5p and (j7.6p . 

/ (^p|u|2 + ip|F|2 + ^-(p^-7p + 7-l)) dx 

Jr«\^ 2 7-1 y 

+ / / (/x|Vu|2 + (/i + A)|divup)dx(is (7.8) 
Jo Jr^ 

= j^^ QpoluoP + ^PolFoP + ::^73y'^^° - 7Po + 7 - 1)^ dx. 

From Remark 13. 3|, p{I + E^) : Vu = 0. Hence, from (jl.Scp and ()1.3ap . we have 

dt{ptiE)=0. (7.9) 
Therefore, from (|7.8p . (|7.9p and the conservation of mass (|1.3ap . we finally arrive at 

/ (Llul' + LlEl' + ^ip-r-^p + ^-l^dx 
Jr3 \2 2 7- 1 J 

+ [ [ (^|Vu|2 + (/i + A)|divu|2)dx(is (7.10) 
Jo Jm.3 

= [ ( ^Poluol^ + ^Pol^oP + -^(pI -7Po + 7-l))dx< R^. 
Jm^ V2 2 7-1 J 

Since /i > is a constant and p G [^) |]i then inequalities (17.ip - (l7.3p follow from (I7.10p . 
and inequality (j7.4p follows from (j7.10p and the following straightforward inequalities: for 
some r/ > 0, we have 

7 1 ( 1^ ^ /r/|x-l|2, if 7>2, 
a; 1 — 7(3; — 1) > < , ,n 

[r/|x-l|2, if < 2 and 1< 7 < 2. 

The proof of Lemma l7. II is complete. □ 

Based on the uniform estimates from Section [6l we have 

Lemma 7.2. Under the same assumptions as Theorem \3.S\ 

ll'5(u) - 1||l°°(0,T;W/i.9) < V^, |k||L°°(0,T;L9) < (7-11) 

for any T G [0 

) ^max] ■ 

Proof. According to (j6.24p . it is obvious to see that 



max ||(T||r9(t) < ^^R. 
te[o,T] 



Hence, we are only left to show 



max ||5(u) — l||wi,9(i) < ^/R- 

tG[0,T] 
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Indeed, for any t G (0, T), we have, by using (jl.3a|) and (|6.19|) . 
I|5(u)(t)-l||2, 

d 



\Po-l\\h + 1^ ^||5(u)(s)-l||2,ds 

|Po-l||2, + ^ r(||5(u)(.)-l 
Q Jo ^ 



a—q 
Li 



X / |5(u)(s) - l|'?-2(5(u)(s) - l)a,5(u)(s)dx)ds ^^'^'^^ 

< IIpo - 1|I2. + - / l|5(u)(s) - l||27'l|9.5(u)||L.ds 
1 Jo 

<S'^ + -i^J^ ||5(u)(s)-l||J-« dsj ^ \\dsS{u)\\l,dsj , 



where 



(5g-6)(p-l) , , 

a = h 1. 

3g-6 



From p.3ap . we obtain 

\\dtP\\LP{0,T;Li{R'-i)) = II^P ' u||Lp(o,r;L9(R3)) + ||/>divu||^p(o,T;L9(R3)) 

< 2||(T||^oo(o,T;L'J{R3))||u||ip(o,T;L°°{M3)) + 2||divu|| j^p(o,T;L9(M3)) (7.13) 

< CR^ + R<CR. 

On the other hand, from the Gaghardo-Nirenberg inequahty, we have 

\\p- ^LiiM?) < c\\p- i|Il2(ir3)||v(/)- i)||^;'(^ig3) < c\\p- i|Il2(ir3)||o-||^7^ 



with e = gl^. Thus, by Holder's inequahty, (fTa . and (|6TT9]) . one has 



which, together (fTT2]) and ([713]) . yields 

||5(u)(t) - < Ci?. 
Hence, according to (j6.19p . we obtain, by letting R be sufficiently small, 

niax max{||5(t;)(t) - l||vyi,9(R3), ||<T(t)||L<,(R3)} < \/^. (7.14) 

te[o,T] ^ ' \ 'J 

The proof of Lemma 17.21 is complete. □ 
Lemma 7.3. For each 1 < / < 3, ^ satisfies 



a*— + u • V— = -— -VE + V { — ] E + Vu— + V— 7.15 

oxi oxi oxi \ axi J oxi oxi 
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in the sense of distributions, that is, for all ip G Cq°{Qt), we have 

/ —dtipdxdt+ / / div(uV')T^ 

-_r [ + V f— ^ E + Vu— + V— ^ 

Jo JR-i V 9xi \dxij dxi dxj 

for any T G (0 



ipdxdt, 



Proof. The proof is a direct application of the regularization. Indeed, one easily obtains, 
using ()1.3cp . 



dtiSsE) + u • V{SeE) = S,{dtE + u-VE) + u- V{SeE) - Seiu ■ VE) 
= SeiVuE + Vu) + u • V{SeE) - Se{u • V^). 

Differentiate ()7.16p with respect to x/, we get 
dt + u • V ' 



(7.16) 



dxi J \ dxi 

Notice that 

\dxij V ^^«/ 
According to Lemma 6.7 in ^24j (cf. Lemma 2.3 in [19]), we know that 

oxi \dxi J \oxiJ V 9xiJ 

in Li(0, T; Li{R^)) as e ^ 0. Hence, 

^-(u-ViSeE)- Se{u-VE)^ ^0 

in L^{0,T;Li(R^)). Thus, letting e ^ in ([TTTll . we deduce 



(7.17) 



dt— + u • V— = -— -VE + V ( — ) E + Vu— + V — , 

OXl OXl OXl \ OXi J OXl OXl 

in the sense of weak solutions. The proof of Lemma 17.31 is complete. □ 
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Using (j7.15p . formally we have, 

dE 



I 



(-) 


dE 


\dxij 


dxi 



dx 



■dx 



I ( 



I 

dE du f du\ dE du 

u ■ V— - — -VE + V { — ] E + Vu— + V — 

dxi dxi \ dxi I dxi dxi 



dE 



dxi 



q~2 



< C[ IIVull^oo ||Vi?|| f „^^3^ + ||ii^||J;^oo(Q^)||u||yl/2,q||V£'||^q^.Jg3^ + ||u||p^2,q ||VS 



< C( llVul 



|Vu||^,.^3, + \/^||u| 



dE 
dxi 

19-1 



dx 



^2,q ||V£'||^,^^jg3^ + ||u||^y2,q IIV-E 



(7.18) 

We remark that the rigorous argument for the above estimate involves a tedious regular- 
ization procedure as in DiPerna-Lions [8] , thus we omit the details and refer the reader to 
[8]. Using (f7J8]l . one obtains 

p 



dE ^ , 



LI 



dE{0) 



dxi 
dE{0) 



LI 



/■* d 




Jo ds 


— is) 

dxi 



ds 



LI 



< 



dxi 
dE{0) p 



p 

+ - 

LI 1 Jo 





dE 




dE 




dxi 


LI Jm.^ 


dxi 



q-2 



dxi 



+ \/^||u| 



Li 





dE 


P-Q |- 


/' 

/o 


dxi 


Li 



iVul 



\VE\ 



IL9 



+ ||u||t^2,9||Vii^ 



IL9 



dE 
dxi 



WE 



ds 



5, 



\Q 
\Li 



dE 
dxi 



(s) I dx 



ds 



< 



dE{0) 



dxi 



Li 



lv^l|^;'[||Vu|u.o||v^iii, 



<6'^P + C 



ds 







+ (1 + VR)\\u\\yy2,g]ds 



max IIV^II + VR+ 1 

te[o,T] 



^ ^ RP I max II VE^II + \/^ + 1 ] . 

t€[0,T] 



(7.19) 



Taking the summation over / in (j7.19p and taking the maximum over the time t, one has, 

- ) RP ( max II V£;|| +VR+i) , 



max ||V^||P <6^P + C 
te[o,T] 



and hence, by letting R, 6 be sufficiently small and using (I6.25p . we obtain, 

max llVSf < + CRP < (/R)^. (7.20) 

te[o,T] 

We are now left to deal with the quantity ||-E||2,9(ir3). To this end, from the Gagliardo- 
Nirenberg inequality, we have 

11-9 



E 



LI 



^ c'II-^IIl2(r3)I|V-E||^, 
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with 9 = -F^rg- Thus, by Holder's inequahty, (j7.3p . and (|6.25p 



l-Ell (5jr-6)p 

L 3<7-6 {0,T;L9 



I ^ -^lll,P(0,T;L9 



lL°°(0,T;L2 

Hence, we have the following estimate: 

Lemma 7.4. Under the same assumptions as Theorem \3.Si it holds 

I|-E^IIl°°(0,T;L'J(R3)) < ^/R, 

for any T £ [0,rmax]- 

Proof. By (fO^ . (lOHl) and ([720]), and letting 

{5q-6){p-l 



< CR. 



a 



3q-6 



+ 1, 



one obtains, 

wmwh 
= \\Emh + 

= \\m\\h + 

= \\Emh + 

< 11^(0)112.+ 







a 

Q Jo 

a 

Q Jo 

a 

Q Jo 



\E{s)\\l-'^ 

a-l 



\E{s)\'i-^E{s)dsE{s)dx] ds 



\E{s)r'E{s) 



VuE + Vu - u • VE 



(7.21) 



(7.22) 



dx ds 



\E(s) 



2\\Vvl\\l^\\E\\l'} + llVul 



Li 



ds 



<\\Em%+-^( f \\E{s)U''^ 



dt 



£-1 

V 



(7.23) 



|u||lp(0,T;VF2.' 



X 2 sup ll-E 

V te(0,Tmax) 



<\\Em%^y—)R[j^ \\E{s)W" dt 

<\\Emh+cR\\E\r,,\_,,^ 

L ii-J-H (0,T;L9 

Then, according to (|7.21|) . one has, for all t £ [0,T„ 



\E 



LI 



< 6^" + CR'' < Vr" 



LI 



+ 1 



(7.24) 



if R is sufficiently small. Thus, (|7.22p follows from (]7.24p . The proof of Lemma 17.41 is 
complete. □ 



Lemma 17.41 together with ()7.20|) and Lemma 17.21 gives 

max max{||5(u) — l||^yi,9(t), ||cr||iq(t), ||T(u)||^i,9 (t)| < CR < 

te[o,T] 



(7.25) 
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Similarly, we can obtain 

max max{||5(u) - l||^yi,2(t), ||fT||L2(t), ||r(u)||v^^i.2 (t)} <CR< Vr. (7.26) 

7.2. Refined estimates on Vp and VE. In order to prove Theorem 13.21 we need some 
refined estimates on || Vp||l2(o,T;L9(r3)) and ||V£'||i2(o,T;L'3(R3))- 

Lemma 7.5. 

for any T £ {0,T^a.x). 

Proof. Taking the divergence in (ll.3bp . one obtains 

AP{p) = diY{diY {pEE'^)) + div(div(p£;)) - div(pu • Vu) - dw{pdtu) + (A + 2//)Adivu. 

(7.28) 

Since, div(/>(/ + E)^) = 0, we get 

div(div(,£)) = ^^(,E,) = = -Ap 

in the sense of distributions. Hence, (j7.28p becomes 

AP(p) + Ap = divdiv(/)SS^) - div(pu • Vu) - div{pdtu) + (A + 2^)Adivu. (7.29) 
Hence, one obtains, using L"? theory of elliptic equations and Taylor's formula, 

I|V/?||l2(o,T;L'J) < C(^||puVu||i2(o,T;L9) + \\dw {pE E'^ )\\ ^2 ^q^t.^,) 
+ IIMu|Il2(0,T;L9) + l|Vu||i2(o,T;L'J)^ 

< C'(ll/'llL°°(QT)ll^"lli2(0,T;L°°)l|u||L°°(0,T;L'J) + W"^ PWl^ {0,T;L9)\\E\\ic^ (^q^) 
+ (Qt) 11^-^11^2 (0,T;L9)I|-S||l°°{Qt) 

+ \\p\\L°°{QT)\\^t^\\L^{0,T;Li) + II Vu||l2(0,T;L'3) ^ 

< C'(^||p||ioo(Q^)||u||i2(o^T;iy2,g)||u||ioo(o,T;L'?) 
+ II V/?||l2(o,T;L9) ll-^lli°°{0,T;VKi.'?) 

+ \\p\\l°°{Qt)\\'^^\\l^{0,T;Li)\\E\\l^(o^T;W^'1)) + -^J 

< C7(i?2 + i?) < Ci?. 

The proof of Lemma 17.51 is complete. □ 

Ir order to refine ||Vi?||j;^2(o,T;L9{R3)), we need the following estimate: 
Lemma 7.6. 

l|5tU||i2(o,T;L9(M3)) < C-RV, (7.30) 
for any T G (0,Tjnax) o-^d some 6 £ (^) !]• 
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Proof. We first notice that, by the Gagliardo-Nirenberg inequality, for q £ (3,6], 

I|u||l2(o,T;L'?(R3)) < C'||Vu||^2(o,T;L2(]R3))||u||^2[|],T;L2(K3)) < Ci?^+^, (7.31) 



with 6 = ^^^^ G (^) Next, we multiply ()1.3bp by dtu and integrate over x (0, t) to 
deduce 



dxds 



JIR» 



JIR» 



div (pE') udxds 



p|5iupdxds+ / (;u|Vu(t)p + (A + /i)|divu(t)p)dx 
= / (//|Vuop + (A + /i)|divuo|^)dx- / / pu • Vu • 9iU 
+ / / div(p£;£;^)atudxds + 

:= / (/i|Vuo|2 + (A + /i)|divuo|^)a!x + Vl, 



with the following estimates on Ij (i = 1...4): recalling Qt = x (0,T), 

1^1 1 < ll\/p9fU||2,2(Qy)||^||ioc(Qj,)||u||^oo(o^T;L2(R3))||Vu||i2(o^T;L° 

< - / / p|atup(i2;ds + Ci?^ 
8 Jo Jm3 



V Pdfudxds 



(7.32) 



o Jo Jr3 

f I p\dM^dxds + CR^- 

o Jo JM3 

Isl < C||Vp||i2(Q^)||£;|||^(Q^)||^9tu||i2(Q^) + (:7||E||i^(Q^)||V^||i2(Q^)||^9iu||i2(Q^) 

< - / / /o|5tu|2dxds + Ci?^ 
8 Jo Jr3 



dt{pE)Vvidxds 



+ 



+ 



p{T)E{T)Vu{T)dx 



poEoVuodx 

< i\\p\\L--{QT)\\9tE\\L2^Q^) + ||^||L-(Qr)ll^t/'llL2(QT))ll^u|lL2(QT) 

+ Ci?3 + (||vp(r)||^2(i,3)||i?(r)||^oc(K3) + ||vi?(r)||^2(K3)||p(r)||^.o(K3))||u(r)y^ 

where, for the estimate I4, we used equations p.3ap . (|1.3cp . Lemma 17.11 and estimate 
()7.3ip . Thus, from (j7.32p . one obtains 



aoad 



WdMll^Qr) ^ CR, 
|Vu||ioo(o,T;L2(K3)) < CR. 



(7.33) 
(7.34) 
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Now, we differentiate ()1.3bp witli respect to t, multiply the resulting equation by dtu, and 
integrate it over to obtain 

~ [ p\dtu\'^dx+ [ {p\Vdtu\'^ + {X + p)\divdtu\'^)dx 

^ dt Jk3 Jlg3 

= - 1 dtp\dtu\'^dx — I dtpu ■ Vu ■ dtudx 

— / pdtu ■ Vu • dfudx — I pu- Vdtu ■ dtudx — / VdfPdtudx /„ oc^, 

Jrs Jk3 Jt^3 [l.6i)) 

dt{pEE^)VdtVidx- [ dtipE)Vdtudx 



7 



where using (|7.34|) . we can control Jj (j = 1...7) as follows: 



l^il 



V p-a\dtu\'^ dx 



< 



dtu\\l4Vp\\L4uUs <CR^Vdtu\\l,; 



and 



IJ2I < ||5iu||L6||Vpb3||u||i6||Vu||i6 < R\\Vdtu\\L4Vu\\l4AuU2 

iJal < \\p\\L^\\Vu\\^s\\dtu\\l, < ||p||z.^||u|||,||Vu|||,||Vatu||i. < CR^\Vdtu\\l,; 

\Ja\ < \\p\\L^ML-4Vdtu\\L2\\dtu\\Le < CR\\Vdtu\\l,; 
IJ5I < C\\dtp\\L4^dtnh2 < C\\Vpu\\L2\\Vdtuh2 

< C7||Vp||i3||u||^6||V5iU||i2 < C\\Vp\\Ls\\VuU2\\VdtVLU2 

<^\\Vdtu\\l,+CR4Vu\\l,; 

\M < \\dtph2\\E\\l^\\VdtuU2 + \\p\\L^\\E\\L^\\dtEU2\\VdtuU2 

< R^Vphs ||u||i6 \\Vdtu\\L2 + R\\VuE -U-VE + Vu\\l2 \\Vdtu\\L2 
<R4Vphs\\VuU2\\VdtuU2 

+ Ri\\Vu\\L4E\\L^ + \\VEU4u\\l. + \\Vu\\l2) \\Vdtn\\L2 

< R4Vn\\l, + ^||V9tu||i. + R\r4Vu\\1, + ||Vu||i.); 

\J7\ < \\p\\L^\\dtEh2\\Vdtuh2 + \\E\\L^\\dtph2\\Vdtuh2 

< i^\mn\\l,+R4Vu\\l, + \\Vu\\l,+R4Vpu\\l, 

< ^\mu\\l, + i?2||vu||i3 + iivuiii. + r4vu\\1,. 



We remark that in the above estimates, we used several times the interpolation inequality: 

|VK2.3(R3) < 11/11^2, 2(]J3) 11/11 Jy2^,(R3) 
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for some 9 £ (0, 1). These estimates and (j7.35p imply that, for R sufficiently small, 

~[ p\dtu\^dx + i^ [ \Vdtu\^dx<R^\\Au\\l,+C\\Vu\\l2+R^\\Vu\\ls. (7.36) 
Integrating ()7.36p over (0,r), we obtain that, using ()7.ip . 



ivau 



< CR2. 



(7.37) 



Here we used the estimate 

\\podtu{0)\\L2 < C{\\uo ■ VU0IU2 + IIAU0IU2 + ||Vpo||l2 + IIV^oIIl^) < 5' 
by letting i = in ()1.3bp . Thus, by (j7.33p . (j7.37p and the Gagliardo-Nirenberg inequality, 



we obtain 



□ 



(7.38) 



\\dt'il\\L'2{0,T;Li 

for some 9 G (^,1]- The proof of Lemma 17.61 is complete. 

With (j7.30p in hand, we can now get the estimate for ||Vi?||j^2(o,T;L9 

Lemma 7.7. For the same 9 G (^, 1] as in Lemma 7.6, it holds 

3-9 

I|V-^IIl2(o,T;L9(M3)) < CR 2 , 

for any T G (0 

Proof. Substituting the following two facts 

dtdWE = div(-u • VE + VuE) + Au, dw{pE) = dw{{p - l)E) + divE^, 

into ()1.3bp . multiplying the resulting equation by |divii^|'^~^div£' and integrating it over 
M^, we can obtain 

p d 



< 



+ 



+ 



+ 



divE\\% + ||divi?||i, 
pdtu\divE\'i-'^divEdx + / puVu\dwE\''-'^divEdx 

VP|div£^|''-2div£;dx + / dw{pEE^)\divE\'^-'^dwEdx 
/ div((p- l)^)|div£;|^"2div£;dx 
I div(Vu£;-u- V^)|div£;|''-2divEdx 



(7.39) 



m=l 



where 



Ml < \\p\\L-^\\dtu\\Li\\divE\\%^ ; 



M2 < ||/9||l=o ||u||j;,g II Vu||ioo ||div £'111^5 < i?||u|| yi/a,? 

M3<c||V/>||L,||div£;||i;^ 



|divi?||i;^ 
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M4 < \\Vp\\L4E\\loo\\dwE\\l,^ + \\p\\Loo\\E\\Loo\\VE\\Li\\dWE\\l-'^ 

< [R^VpU. + R\\VE\\L.)\\dwE\\l-'; 

M5 < Hp- l||Loo||V^||L,||divE|||;^ + ||V/)||L,||^||Loo||div^|||;^ 

< R{\\VE\\l. + WVpU^mwEWl','; 

Me < (||Vu||l.c||V£;||l9 + || Au||l, ||£;||L.c)||div£;||^;^ < R\\u\\w2,,\\dwE\\l~^ . 

With those estimates in hand, we multiply (j7.39p by |divi?|^~^ to deduce that, using 
Young's inequality, 

l^lldivi^lli, + lldivSlli, < C\\dtu\\% + R^uWl,,,, + llVpIlL + R^\\VE\\l,. (7.40) 

On the other hand, we still have 

llcui-l^lli, < ||^||ioo||V-E|||, < CR^\\VE\\l,. 



Hence, substituting this into (|7.40p . we get 

|^||divi?||i, + ||Vi?||i, < C\\dtu\\h + R^uWl,,,, + \\Vp\\l,+CR^\\VE\\l,. (7.41) 
Integrating ()7.4ip over (0,T) and using the estimates ()7.27p . ()7.30p . we obtain 

3-9 

I|V-^IIl2(o,T;L<?(M3)) < CR 2 . 

The proof of Lemma 17.71 is complete. □ 
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